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Abstract 

Analytic solution is given in the steady state limit t ^ oo for the system of Master 
equations describing a random walk on one-dimensional periodic lattices with arbi- 
trary hopping rates containing one mobile, directional impurity (defect bond). Due 
to the defect, translational invariance is broken, even if all other rates are identical. 
The structure of Master equations lead naturally to the introduction of a new en- 
tity, associated with the walker-impurity pair which we call the quasi-walker. The 
velocities and diffusion constants for both the random walker and impurity are given, 
being simply related to that of the quasi-particle through physically meaningful equa- 
tions. Applications in driven diffusive systems are shown, and connections with the 
Duke-Rubinstein reptation models for gel electrophoresis are discussed. 
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1 Introduction 



In equilibrium statistical mechanics the number of exactly (analyticaly) solvable models 
is very limited (for a review see [jjl). In nonequilibrium statistical mechanics this number 
is considerably smaller, and no review has been published, to our best knowledge, so far. 
Part of this difference lies in the starting points. In equilibrium physics, one usually begins 
with the well-known Boltzmann distribution. By contrast, arriving at this step is already 
non-trivial, even for steady states, since one must solve, say, a Master equation first. 
To find an analytical solution for such equations, for typical nonequilibrium systems, is 
usually impossible because of the very large number of different equations involved, since 
the unknown is a function of the configurations of the system. However, we can hope to 
find analytic solutions if the dimensionality of our system is low enough. As an example, 
consider the random walk of a single particle on a one dimensional (ID) lattice, with time- 
independent hopping rates Wij, where Wij represents the probability for the particle to 
jump from site j to site i in one step (or per unit time). These rates may be arbitrary, so 
that they are not required to satisfy the condition for detailed balance, so as to describe 
a general, nonequilibrium system. With the additional assumptions (a) that only jumps 
to the nearest sites are allowed {\i — j \= 1) and (b) that periodic boundary condition 
(PBC) is imposed, this model was analyticaly solved in the steady state limit by Derrida in 
1983 (see [^). Explicit expressions for the steady state distribution, velocity and diffusion 
constant were given. However, if, say, jumps to next nearest neighbours are allowed, an 
analytical solution seems to be impossible to give On the other hand, since 1983, 
several other types of ID systems (e.g., @, @) have been solved. 

In the next two subsections of the Introduction we explicitely present on a "checker- 
board" the one-dimensional periodic hopping model without and with a mobile directional 
impurity. Then we give a short description of the main body of the paper. 

1.1 One-dimensional hopping model without impurities 

The hopping model presented by Derrida in |^ can best be visualized as a (A^ + 1) x 1 
checkerboard, i.e., a string of 1 squares). This is filled with pieces (particles) num- 
bered from 1 to A^ (from left to right) and one hole, which will play the role of the random 
walker (Fig.l) . The periodic boundary condition is equivalent to identifying the squares 
at the two ends as nearest neighbours. By saying "that the hole is at site n" we mean 
that the empty square is found in between particles with labels n and n + 1. Since the 
dynamics will be restricted to particle-hole exchange only, we may regard the hole as the 
random walker, and will use the terms'hole' and 'walker' interchangibly. The proposition 
"the hole jumps from site n to site n + 1 in the next time step" on Fig.l is equivalent 
to "shifting the hole from the square in between particles n and n + 1 to the square in 
between particles n + 1 and n + 2" . The rules of evolution are as follows: 

1) a square can be at most occupied by a single particle, 

2) a particle can interchange squares only with the empty square (hole) and only if that 
is a nearest neighbour. 
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Figure 1: The translationally invariant one dimensional hopping model with periodic 
boundary conditions, on a checkerboard. The particles are represented by numbered 
circles (pieces). 



3) the particles are not allowed to interchange positions among themselves, 

4) whenever the hole at site n jumps to the left (to site n — 1), it does with probability 
Wn-i,n, and, for a jump to the right (to site n + 1), the probability involved is Wn+i,n- 

Note that the system above is translationally invariant, in the sense that the jump 
rate of any particular particle-hole pair is independent of its location on the checkerboard. 
In particular, after the hole jumps through the entire system a few times, (effectuating 
complete revolutions) the string of particles will be displaced, relative to the "board" . Nev- 
ertheless, the physics of the evolution remains unaltered. Next, we break this invariance 
by introducing a "kink" in the checkerboard. 

1.2 The model with one mobile directional impurity 

In this paper we introduce a mobile, directional impurity into the periodic hopping model 
above. In terms of the checkerboard, this is a specific "bond" between two adjacent 
squares, such that rates of particle- hole exchanges across it are different (in the sense that 
they differ from the W rates) . Thus, we also refer to it as a "defect bond" (DB). Since 
the string of particles can be displaced relative to the checkerboard, this DB can move, 
relative to the string. In this sense, we use the term "mobile directional impurity" . A more 
transparent way to display the model is represented in Fig. 2, where the DB is shown as a 
kink in the checkerboard. Only across this kink are the probabilities of the particle-hole 
exchange different from the Ws specified above. The motivation for such a kink will be 
discussed below. Here, we simply state the model. In addition to the rules l)-4) from the 
previous subsection, we have one more, specifying the jump rates through the kink: 

5) the probabilities for the interchange of the hole with a particle through the DB (thick 
dashed line) are q, if the hole moves upward, and q', if downward, independent of the 
particle involved. This is a simplifying assumption, allowing less involved formulae. Our 
methods will lead to an analytic solution even if these probabilities do depend on the 
particle. 

Note that, in both models, the order of the particles is never changed, i.e., there is 
no mixing among the particles. Therefore, when the hole is in between particles n and 
n+l, the hopping probability to the left (right) is always Wn-i,n {Wn+i,n), except when 
the DB is involved. However, the string of particles as a whole shifts with respect to the 
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Figure 2: Translational invariance is broken by introducing a kink in the checkerboard 
(lattice). The jumps across the kink (dashed hne) obey different rates, i.e., q and (7', inde- 
pendently of the particle involved. If the particles carry charges, then in the case of an 
external electric field E the jumps across the kink are not biased. 

checkerboard (and the DB), as the hole wanders. 

Fixing our reference frame to the string of particles, which turns out to be mathematically 
convenient, we see that the DB 'moves', if and only if the hole passes through it. Further, 
the jump rates from site n to n it 1 are not always Wn±i,n, but take on the value q or q' 
whenever the jump is across the DB. Naively, one may expect that an analytic solution 
is impossible, given the seemingly time-dependent rates. Nevertheless, we are able to 
find a solution. The key observation is the existence of a specific linear combination of 
the positions of the hole and the DB, with which the system of Master equations can 
be rewritten in a form describing the hopping of a single random walker (called quasi- 
walker) on a chain of length N(N + 1). Exploiting the results of Derrida, the probability 
distribution, velocity and diff'usion constants of the quasi-walker, as well as those for the 
hole and the DB, can be explicitly computed. 

An immediate application for the model with directional impurity is a generalized 
asymmetric exclusion process (ASEP) Q, Q. Consider the pieces as charged particles. 
More precisely, let the particles be charged with Qn and let us apply an external electric 
field E pointing along the chain (See Fig. 2). Neglecting the interaction among charges, the 
hopping rates are solely determined by the infiuence of E on the charges. Away from the 
kink, these rates would be biased, due to E. However, across the kink (or DB), particle 
hops are perpendicular to the external field, so that the corresponding probabilities {q and 
q') are not biased by E. 

Another application, which is also the physical motivation of our model, comes from the 
connection to the reptation model of Duke- Rubinstein for gel electrophoresis Q. In Ref. 
Q we introduced a model for electrophoresis of polymers with impurities, i.e., polymers 
having segments (reptons) insensitive to an external electric field E. We showed that the 
eff'ect of a single impurity is equivalent to having a single kink in the ASEP above. The 
drift velocity of the polymer chain, being related in this model to that of our random 
walker, is therefore directly aff'ected. The details of this application, which is a particular 
case of our model, are presented in the last section of this paper. 

The structure of the paper is as follows: in Sect. 2 we provide a precise formulation for 
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the model with one mobile directional impurity, at the level of Master equations. Here, we 
show that the structure of these equations leads us to a natural definition of a new entity, 
namely, the quasi-walker, which stands for the hole-DB couple. A full solution for the 
steady state probability distribution is then given. From that, the expressions for the ve- 
locity and diffusion constant (in the steady state) of this quasi-walker are derived. Section 
3 is devoted to a detailed description of the original system in terms of the quasi-walker. 
Using a technique from the theory of random walks on rings , see [||, (and referred as 
replication method by us), the expressions for the velocity and diffusion constant for the 
hole and the DB (with respect to the chain) are obtained. The remaining sections deal 
with special cases and applications. 

2 Definition of the model and the quasi-walker. 

In this section, we will present another representation of our model which is equivalent to, 
but simpler than, that shown in Fig. 2. Though there are two moving objects (the hole 
and the defect bond), the motion of the DB is completely determined by the dynamics of 
the hole, so that there is, effectively, only one degree of freedom. This single co-ordinate 
is then associated with the quasi-walker. 

Instead of a chain on a kinked checkerboard, our system can be redrawn as in Fig. 3a)- 
d). Here, we have our fixed reference frame to be that of the particle chain, arbitrarily 
choosing one particle as the "first" in the labeling 1,2,..., N. Representing the particles 
(hole) by solid (open) circles, we see that the hole is located in between the particles, 
i.e., at positions denoted by short vertical lines. The DB (which is the kink in Fig. 2) can 
also be thought of a "wanderer" amongst the particles. We denote its position by a long 
vertical line. The only complication arises when the hole comes "in contact" with the DB, 
i.e., when both of them are located between two successive particles. Then, there are two 
distinct states of the system, depending on which side of the DB we find the hole. As a 
result, we are forced to draw two short vertical lines on either side of the DB. 

To keep track of these positions we introduce the following variables: 

• Associate the position of the DB with b, i.e., when the defect bond is between particles 
b and b + 1. Thus, b G [l,iV]. 

• Let m denote the position of the hole relative to the DB. Since the hole occupies a square 
in the checkerboard, it is clear that this square can be on either side of the DB (shown by 
solid line in Fig. 2). Thus, m assumes one more value than b. We choose m € and 
let m = and N be associated with the positions next to the DB (See Fig. 3). 

For clarity, let us illustrate the dynamics of the hole near the DB. There are 4 representa- 
tions of the same piece of chain at consecutive time steps. In the first (Fig. 3a), the hole is 
at m = — 1 and the bond is at b. In the next step, the hole interchanges with particle 
b, arriving at m = A (Fig. 3b). Suppose it next interchanges with particle 6-1-1. In the 
checkerboard, this would be an exchange "across the DB", so that the DB will now be 
found between particles b + 1 and 6-1-2. At the same time, the hole will be "on the other 
side of the DB". Thus, both the hole and the DB arrive at new positions (Fig. 3c). Finally, 
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Figure 3: An equivalent representation of the model from Fig. 2. The solid circles represent 
particles. The short vertical lines mark the possible positions of the hole represented by 
an open circle. The DB is shown as a long, solid, vertical line. As the hole passes through 
the DB, the latter encounters a one-step shift in the same direction, a)-d). 

according to our scheme, we must "relabel" m accordingly to the value m = 0. In the 
last frame, the hole moves further to the right. The DB is not affected, so that only m 
increases to 1 (Fig.Sd). Of course, if the hole moves to the left, then we must reverse the 
previous process. Let us emphasize that, in the single move between Figs. 3b and 3c, b 
increments by 1 while m drops from A'^ to 0. These rules can easily be verified by playing 
on the checkerboard. 

2.1 The system of Master equations 

The full probabilistic information about our system is contained in P{b,m;t), the proba- 
bility to find the DB and the hole at b and m, respectively, at time t. The time evolution 
of P will be controlled by the jump rates of the hole, defined in rules l)-5) above: Wi_j, g, 
and q'. Recall that the subscripts on W refer to the positions of the hole being between 
particles i and i + 1, etc. In terms of our variables b and m, this position is clearly just 
b + m, modulo A^, since our chain is a closed one of N particles. Thus, it is natural to 
introduce the notation: 

^m,m' = Wb+rn,b+m' , (1) 
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with the understanding: 

Wi+N,j+N = W,,j . (2) 

Since the hole can only exchange with its neighboring particles, the restriction \ m — m'\= 1 
applies in (|^). Combining these remarks and following the special rules associated with 
the DB (Fig. 3), the evolution of P is given in terms of a system of Master equations: 



m = 0: dtP{h, 0) = P{b - 1, N)q + P{b, 1) Wq^ 

-P(6,0)(g' + W^y 
m / 0, iV : dtP{h, m) = P{b, m - + P{h, m + 

m = N: dtP{b, N) = P{b, N - 1)W^n,n-i + Pib + 1, 0)g' 

-P(6,iV)(<_i,^ + g 



(3) 



where we have suppressed t for simplicity. 
2.2 The steady state limit 

One of our central interests is the steady state limit. Assuming that it exists, let us define 

= lim P{b, m; t) . (4) 



t— >oo 



Then, from (0), these satisfy the following set of linear equations: 



m = N: + X',+\' -X%(^q + = 



These equations can be solved by introducing the steady state probability current densities: 



= X^ W'' , - X^ ,W'' 

— ^'-m" m—l,m ^'-m—l"m,rn—l 



(6) 



Now, one can rewrite (Mj as: 



m = 
m = N 



X 



6-1 

N '-- 



'-■m+1 1 
b _ vb+lJ n 



(7) 
(8) 
(9) 
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From (|8|), we find 

Cat = c^r-i = ••• = C2 = c5 . (10) 

Examining the right hand sides of (^) and we can go one step further: 

4 = cr . (11) 

In other words, none of these c's depend on b or m, so that we might as well write 

cl = c . (12) 

This equation can be easily interpreted as: the probability current density in the steady 
state is a constant. This result is typical for systems with a periodic, one dimensional 
configuration space, where the current has only a single component and can be fixed by 
the conservation law. Since we have two variables (b, m), we may regard our configuration 



space as "two dimensional", so that (12) is perhaps surprising. The resolution is best dis- 
played by explicitly writing the entire system (|^, |8|, |^ using (|^), as in Fig. 4. Note that, 
between the lowest and highest equations, we have the periodic boundary condition (^), 
i.e., = Xj^, and Xq~^^ = Xq. Thus, the whole system of equations can be regarded 
as those for a single, periodic chain. Explicitly, we introduce: 

Xk=Xi, be[l,N], mG[0,N] , (13) 

where 

k = l + {b-l){N + l) + m (14) 

simply numbers the equations of Fig. 4 starting from the bottom. Thus, k £ [1, M], with 

M = N{N + 1) . (15) 

Observe that k uniquely determines the pair (6, m) and vice- versa: 

m = (fc-1) mod(A^-M) , (17) 
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where [x] represents the integer part of the real number x. To display the structure of this 
single periodic chain more clearly, let us introduce: 

Wk,k-i^w!^^^_^ and Wk-i,k^w!^_,^^ , (18) 

for m £ [1, N]. For m = 0, we first define 

kb = l + {b-l)iN + 1) . (19) 

Referring to (^, ^), we see that 

Wk,,kt-i = q and = q' , (20) 

for all b. Finally, to stress the M-periodicity of these rates, we write 

Wk+M,l+M = Wk,l . (21) 

With this notation the system from Fig. 4 reads: 

XkWk-i,k - Xk-iWk,k-i = c , kG[l,M] . (22) 

Recast in the form of (p^, our system is exactly the one studied by Derrida, |p. The 
only difference is that, in our case, the Xk stands for the steady state of the "defect bond- 
hole" pair, as opposed to a single particle as in We will refer to this new entity as a 
"quasi- walker " , and treat it as a single random walker on a chain of length M. 

Lastly, we remark that not only the steady state limit can be reduced to the case in 
Q, but also the full time dependence of our system (^) can be mapped. Defining 

Pk{t) = P{b,m;t) , (23) 

the system of Master equations (^ can be written as: 

dtPk = Pk+iWk,k+i + Pk-iWk,k-i - Pk (Wk+i,k + Wk^i^k) . (24) 

2.3 Application of Derrida's results to the quasi-walker. 

In this paragraph we apply directly the results from Q to obtain the steady state proba- 
bility distribution, the velocity and diffusion constants of our quasi-walker. As expected, 
all quantities are rational functions of the jump rates. Thus, 

Xk = Frfc, (25) 
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Figure 4: The rewritten system. 6 labels the blocks, m the corresponding equations. 
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where 



1 



and r is a normalization constant: 



M-l i 



i=l j=l yyk+j+i,k+j 



(26) 
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k=l 



The velocity and the diffusion constant are given, respectively, by: 
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(30) 



While it is straightforward to arrive at (H-l^), a subtle process is involved in finding 
( p^j29| ). Before proceeding, let us review this method which we call "replication". 
Instead of studying a finite chain with periodic boundary conditions, replicate it into an 
infinite, periodic chain. Thus, we allow, for example, the indices i,j in Wi^j to be any 
integer, together with the periodicity constraint, (^). Similarly, the distribution Pk{t) 
in (^) is now defined for all k G [—00,00], and satisfies (p^) with (]2l|). To compute the 
steady state quantities defined for the original chain, we must sum over all equivalent sites 
on the replicated chain. On the other hand, the velocity and diffusion constants can be 
found from the first and second order moments of Pk- Introducing 
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(31) 



the velocity is naturally 



(32) 
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while the diffusion constant is given by: 

D = \im ^[X2{t) - {x{t)f] . (33) 
t^oo dt 

Note that x{t) is the position of the quasi-walker at time averaged over all paths. In 
Sect. 3 we present in more detail an extension of this replication technique to describe the 
hole and the defect bond separately. 

Before proceeding, let us present these results in terms of the original jump rates Wij, 
q' and q. For convenience, let us first define the following expressions: 

Sib,m)^ Y: n , me[0,N-2] , (34) 

5(6,A^-1) = , (35) 



r(6,m)^f:n^±n^^^^^±I^, mG[l,iV] , (36) 

m—l 



T{b,0)=0 . (37) 



The last term of (|34D and (3£) also appear frequently, so that we define: 



Xs{b,m)^ n 7;+-+'-^'"+-+' , mG[0,iV-2] , (38) 

1=1 >^b+m+l+l,b+m+l 

Xs{b,N -1)=0 , (39) 



\ /L \ TT ^b+m-l,b+m-l+l 

XT{b,m) = [[— 

1=1 '^^ b+m-l+l,b+m-l 



m £ [1,N] 



(40) 



At(6,0) = 1 . (41) 

Observe that the above quantities depend on neither q' nor q, and that the W^s satisfy 
periodicity As a result, they are also periodic: e.g., S{b + N,m) = S{b,m), etc. 
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Finally, let: 



and 

, N 



After some rather tedious algebra, we arrive at the explicit expressions: 



where 



Wfcb+m = T7} |f/[l + m)] + g'AT(6, m)ufc6_i+7v| 

m G [0, iV - 1] 
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The normalization constant (27) now takes the form: 
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>^sib, m) 

Wb+m+l,b+T 



From (|2^), we see that the current c can be related Q to the velocity of the quasi- 



via : 
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'm ' 
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with V given by (|28|). Verifying 

A 

Wk+i,k V g 

we have: 

V = TMU . (51) 



In other words, U from ( p3| ) contains the essential information on the velocity y. Further, 
it indicates whether our system satisfies detailed balance or not, i.e., if we have an equi- 
librium or a non-equilibrium steady state. The former is the case if and only if c = 0, so 
that both U and V vanishes. The condition to stop the motion can be summarized in a 
simple equation: 

q/q' = \ . (52) 



Its interpretation is clear. The left hand side is the bias in the jump rate across the defect 
bond, while the right hand side represents the bias in the opposite direction across the 
entire chain of particles. When these balance, we have equilibrium. 

In the next section, we will see that the expression (^ij) for the velocity of the quasi- 
walker plays a central role in determining the velocities of the DB iybond) and the hole 
{Vhoie) with respect to the chain, as well as the relative velocity between the two {Vrei)- In 
particular, an extension of the method of replicating a finite closed chain will be shown in 
the derivation of these quantities. Following this method we also can derive the diffusion 
constants for the bond and the hole Dhond and Dhole in the same steady state limit. 

3 Replicating the finite, closed chain 

In a finite, closed chain, a hole will visit the same site infinitely many times in the t — >■ oo 
limit. As a result, it is not easy to find asymptotic properties like steady state velocities 
and diffusion constants. To circumvent these difficulties, the method used (see for e.g. |^, 
0) is to replicate the closed chain into an infinite, periodic chain (with periodic jumps 
rates). Since there are infinitely many sites, though many are "equivalent", this infinite 
chain greatly facilitates the computation of asymptotic quantities. Here, we will generalize 
this method and arrive at expressions for the velocities and diffusion constants of both 
the DB and the hole. The extra complication in our case is that, instead of having one 
random walker, we have two moving objects: the hole and the defect bond. Therefore, we 
have to define our replicated model carefully so as not to lose any physical information 
about the finite, closed chain. 
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Figure 5: The replicated, infinite chain. The long, solid vertical line represents the DB's 
position while the long, dashed vertical ones show its images, j labels the number of 
complete blocks (containing A'' particles) away from the DB. 

On the infinite chain, the position of both the hole and the DB can take arbitrary 

values, so that they can be found with any distance in between. Let us define the position 
of the DB, relative to some arbitrarily chosen "origin", to be (3 and the distance from the 
bond to the hole to be /i. Analogous to b and m, these differ in that [— oo, oo] 

Now, the dynamical rules associated with these are not so simple, however, since they 
must reflect the fact that, on the closed finite chain, the DB "advances" by a step each 
time the hole moves forward N steps. To simulate these advances, we introduce an infinite 
set of "images" of the bond, located at multiples of from the DB. These images serve 
as "condition lines" for the motion of the DB itself, i.e., whenever the hole passes through 
one of these images, the fi changes by unity. Of course, in accordance with the idea of 
an image, all the images also move correspondingly. In this way, though the two objects 
can be arbitrarily far apart on the replicated chain, the properties, or the structure, of the 
original chain are replicated faithfully. Referring to Fig. 5, it is clear that there are unique 
integers i and j, such that 

l3 = b + iN and /j, = m + j{N + 1) . (53) 



Note that, as expected, b is just P mod (N). Now, j represents the number of complete 
replicas of the original chain lying between the DB and the hole, while m is just the 
distance, as in the previous section, from the hole to the nearest bond-image to its left. 
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Let us caution the reader on an extra complication associated with the second equation 
in (|53|) . The reason for the factor {N + 1) hes in our having to distinguish between a hole 
immediately to the left/right of all the images of the DB. Therefore fi represents the number 
of states or positions between the bond and the hole. However, the "actual distance" , in 
terms of number of particles (or units of lattice constant), between the DB and the hole 
is only m + jN = fJ- — j (see Fig. 5). In other words, while it is necessary to associate two 
positions with each image for the hole, there is only one "lattice spacing". This subtle 
distinction will be crucial when we compute the velocity and diffusion constants, both of 
which must be based on the physical distance. In particular, the distance of the hole from 
the origin of the infinite chain is simply P + fi — j. 

To distinguish the statistical mechanics of the replicated chain from the finite one, we 
will use notation P{P, /u; t) for the probability distribution. To write down the Master 
equations for P , based on the original set (^), we first transcribe W^^,, the transition 
probability for hole-jumps on the replicated chain. Using the periodicity property (^), we 
see that these are simply given by 

, =W^ , . (54) 

" m,m' V"^y 

With the notation /Xj = j{N + 1), and suppressing t, we define the currents 

m / : J^iP) ^ P(/3,/x)<_i,^ - PiP,l^- l)<^_i, 
^ = 0: J|^(P)^P{P,^,J)q' -P{P-l,fi,-i + N)q . (55) 



In terms of these, the system of Master equations for P{(3, /i) reads as: 



m 







m 



N : 



dtP{P,fij) 
dtP{P,fi) 
dtP{f3,iMj+N) 



JLi{P)- 



(56) 



3.1 Positions, velocities and diffusion constants 

To reach our goal of computing velocities, etc., we need some notion of positions. For 
that, we introduce the following first moments of P: 

oo oo 
A»=-oo /3=— oo 
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oo oo 

{m= E E , (58) 



OC oo 



m)= E E jPif^^i^) ■ (59) 

/i=— oo /3=-oo 



Is easy to see that, the average position for the DB is just {f5{t)). For the average 
distance between the hole and the DB, we define 

{5{t)) = {i,{t))-{j{t)) , (60) 
so that the average position of the hole is 

m) ^ m) + {s{t)) = m)) + im - m) ■ m 

The time derivatives of the moments defined above can be used to compute velocities: 

Vbondit) = J^m)) , (62) 

v,{t) ^ |(/x(t)) , (63) 
= J^m) ■ (64) 

The first clearly represents the average velocity, with respect to the chain, of the DB. 
Naively, the second can be thought of as the velocity of the hole relative to the bond on 
this infinite chain, since /j, plays the role of m here. However, in view of the distinction 
between fi and the "actual distance" ^ — j,we emphasize that the "actual velocity" of the 
hole relative to the bond, is 



v,,i{t) = -^{6it)) = v^,it) - vj{t) , (65) 



while the hole's velocity, relative to the chain itself is 
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Vhole = ^(^(*)) = Vbondit) + Vrelit) ■ (66) 



The diffusion constants for tlie DB and for the hole can be found similarly: 

dbondit) = \[j^iWHt))) - j^mt)f)] , (67) 



dhoUt) = l[j^{{h'm-j^imf)] , (68) 



where 



oo oo 



|[(/5'(i))]= E E (^'^tP{P,^.) , (69) 

Ai=— oo /3=— oo 



oo oo 
At=— oo /3=-oo 



d 



Needless to say, all the six second order moments are needed (all quadratic combinations 
of (3, ^ and j). In this paper, we will show explicit calculations for the first order moments 
and velocities only, since the computations for the second order moments are considerably 
lengthier yet providing no additional insight. 

3.2 Reductions to the finite, closed chain 

As the original system consists of only a finite, closed chain, it is important to make the 
connection between the replicated, infinite chain above and the finite one. The simplest 
way is to note that i and j play the role of replica indices. Thus, we have: 

oo oo 

P(6,m;t)= E E ■ (71) 

i=— oo j=—oo 



To facilitate the discussion of the moments (57-59), we define 



oo oo 



</.(&, m)= E E /3^(A/^) , (72) 



i=— oo j=— oo 
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oo oo 



i=—ooj=—oo 

so that the positions and velocities are given by 

N N 
6=1 m=0 



(73) 



(74) 



TV N 
6=1 m=0 



(75) 



N N 



(jit)) 



N + 



1 



6=1 m=0 



+ 1 



N N 
6=1 m=0 



(76) 



N N 

Vbondit) = E E dt4>ib,m) 

6=1 m=0 



(77) 



and 



TV N 



^(*) = E E ^*^(^' 



m 



6=1 m=0 



(78) 



1 



+ 1 



TV TV 



^w-EE"^w, 



6=1 m=0 



(79) 



To compute ( [77y79|) and the time derivatives of the six second order moments, we should 
use the Master equation for P{b, m) to eliminate dtP in favor of the currents and derive 
similar equations for (f)(b,m) and ip{b,m). Not surprisingly, these equations will comprise 
of current-like terms. Letting a be any of the quantities P, cp oic ip, we define: 



18 



m = 0: J^{a)=a{b,0)q' -a{b-l,N)q 



(80) 



and note that these represent the net probabihty current density from state (6, 0) to 
(6 — 1, A^) for m = and {b,m) to (b,m — 1) for m ^ 0, respectively. Using these, the 



definitions d7l-73) and the Master equation for P, we obtain: 



m = 
< m < iV : 
m = N : 



dtPib,0) = J,\P) - J,\P) 
dtP{h,m)=jt^^{P)-jl^[P) 

dtP{b,N)=j'^+\P)-jUP) 



(81) 



m = : 
< m < iV : 
m = N : 



dtHb, 0) = J^{(t>) - Jo'(0) + P{b - 1, N)q 
5i<A(6,m)= J^+,(0)- J-^(<A) 
dtHb, N) = J^^\<P) - J'k{4>) - P{h + 1, 0)g' 



(82) 



m = : 
< m < iV : 

m = N : 



dMb, 0) = - Jo'(^) + P{b - 1, N)q - P{b, 

dMb,m) = j:^+i(V) - Ji{i^)+P{b,m- 1)h1,„„i - 
P(6,m+l)iy^^„+i 

dMb,N) = Jo'^^(V) - J^(V') + P(6,iV - l)VF^,^-i - 
P(6+l,0)(?' 



(83) 



Note that ([81]) is just (H), as it should be. Now, keeping in mind the periodicity property 



P{b + N,m) = P{b,m), we may express for the velocities (77- 79) in terms of the proba- 
bility current-density alone: 



N 



Vbond{t) 



b=l 

N N 



b=l m=0 
Vj{t) = Vbond{t) . 



(84) 

(85) 
(86) 
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The last of these, obtained by exphcit computation, can be understood intuitively. Since 
j measures of the net number of times the hole traverse around the (finite) chain, while 
the DB moves one step each time the hole crosses it, the displacements in j and the bond 
are locked. Prom these, we easily arrive at 

Vrel{t) = Vfj,{t) - Vj{t) = Vf,{t) - Vbondit) (87) 

and 

Vholeit) = Vf,{t). (88) 

The f-derivative of the second moments, (|6^) and (^) can be computed similarly: 

, N N 

-{{P\t))) = -2Y: J^{<P) + E[^(^' + 0)q']. (89) 
f,=i b=i 



and 



, N N N N 

j^{{h\t))) = -2^ E Jt{ct>) - 2^ E 

6=1 m=0 b=l m=0 

2 N N N 



-EE ^^tiP) + 2 E ^(^ - 1' 



^ 6=1 m=0 6=1 



N N 
+ E { [P{b. 0)q' + P{b - 1, N)q] + E [Pib, vfi)Wl_^^^ + P(6, m - l)t^^,^_i] }. (90) 

6=1 m=l 



Note that these expressions involve sums, as opposed to differences, of terms like Pq and 
PW ^ so that they are not simply related to the currents. 

These results are general, applicable for all t. Our interest lies in the long time steady 
state limit, which is the focus of the next subsection. 

3.3 Steady state limit 

In the t — > oo limit, we already know that P(6, m\ t) X'^ (Section 2.2-3). For the other 
two, our expectations of the asymptotic behavior lead us to the ansatz: 



0(6, m) ^ + , 

^lj{b,m) ^ Xlt + ^l . (91) 
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To determine the four unknowns, we insert these asymptotic forms into (81- ^3|) 
obtain a systems of equations for them. Letting a denote either uj or A, we have 



m = 
N 



m 



For $i and ^t,, we find 









and 



(92) 



m 







m 7^ 0,iV 
m = N 



4 



^0 1 



(93) 



and 



m = : 

m 7^ 0,iV : 
m = N : 



An 



J'' 



Jo'(^) + X^-1q 



Tram 



m— 1 ' ' m.m- 



l^N,N- 



m+1 ' ' m,m+l 

Aq q . 



(94) 



Now, ( |9^ ) is the same set of equations as X^, so that u;^ and Aj^ must be both propoer- 
tional to Xt, ■ Let us write 



\b _ T/ 



(95) 
(96) 



where the constants Vbond and can be fixed by the normaUzation condition 



EN sr^N yb 1 . 
b=l l^m=0 — ^■ 



N N 

Vborid = ^^^m ; 
6=1 m=0 

N N 
6=1 m=0 



(97) 
(98) 



Using (77-^) and (91), we arrive at the meaning of these constants: 



Vbond = lim -jriPit)) = lim Vbond{t) , 



(99) 
(100) 
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In a similar way, we have Vj = \\Tat^ood{{vj(t)))/d't = Vbond-, Vrei 



V, 



Va. 



hole — ' fi 

Next, we apply the results from Sect. 2.2-3, namely: 



lim J^{P) = Ji{X) = c 

t—>QO 



V 

'm 



Va - Vbond and 



(101) 



and find, from (|84| - p5D , 



Vbond = -Nc 



V 



Vrel 



V 



N + 1 ' 



= -Mc = V 



V 



N + 1 
N 



N + 1 

Vhole = V 



-V 



(102) 
(103) 
(104) 

(105) 
(106) 



Note that these results can also be obtained, according to ( |97| ) and (|^)) by summing ( |9^ ) 
and (^ ) over m and b. 

Finally, we turn our attention to the diffusion constants ( |67| - |6^) . Using the linearity 
property of i7m(cK) as function of a , as well as the expressions (|95n96|) and (99 - 100 ) with 



(P4-B5), we verify that terms linear in t cancel. The expressions are somewhat simpler if 
we used the renumbering scheme ( [l^jT9| ) for the quasi-walker. Accordingly, all quantities 
of the form are now denoted by A]^, where k = 1, ..,M. In this notation, the current- 
density type quantities in ( |80|) are summarized in one line: 



Jk{A) = AkWk-i,k - Ak-iWk,k- 



k = 1,..,M 



(107) 



As in the previous section, we can solve these for Ak- If we sum over k, an important 
relation is obtained: 

^ _ / rM\ ^' — - 

Y^Ak=[-—)Y.^kJk+M) , (108) 

k=i \ V / i^^i 

where Uk is given by (^). Using this notation, the diffusion constants are: 
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N ^ N N N 

Dbond = - E ^0 (^) + 2 E + ^n'i) - ^bond E E ^ 

6=1 b=l b=l m=0 



b 

m ! 



(109) 



M M N 

Dhole = - E W - ^ E ^-fc +T.^h-1Q 
k=l k=l 6=1 

M MM 

E w - ^ E *fc + E Xk-iWk,k- 

k=l k=l k=l 

M _ y 

^ -ErnX,--.^^ • 

fe=i 



+ 



N 
N + 1 

V 



(110) 



Now, we need to find expressions for Jk{^) and Jk{^)- These can be found by solving 
(1) and (11: 



y 1 ^ , V m 



jkm 



N + 



+ IM 



M 



N + 



^ TTf E ^'^k' + Sk,ki,Xk^^iq + c$ 
k'=i 



(111) 



and 



__ V ^ — 

Jk{^) = T7 E ^'^fc'+fc-l + Xk_iWk,k-l - T7 E ^'^fc' + 
A;'=l fc'=l 



(112) 



where m = (A; — 1) mod {N + 1) and c$(^) are constants which will disappear from the 
final expressions. Using ( |108D and substituting (111, |112| ) into ( p,09| , |11C| ), we finally arrive 
at explicit forms for the diffusion constants: 



bond 



{N + iy 

TV 



MM M 

^ E E ^'^fc'+fc + ^ E ^fcb-i^fcb-i' 

k=l k'=l k=l 

M 



(iV + 



Uk-1 - rk m 



k=l 



M + 2 V 
2 (iV + l)2 



(113) 



and 



/MM M \ 

Dhole = h/ E ^fc E ^'^k'+k + M E uk-fkWk+i,k - V 

\ k=l k'=l k=l I 



M + 2 
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+ T'^N'^\^Uk^-irk^^iq -Uk-irk^iWk,k-i + - ^(ufc_i - rfc)m. (114) 
k=i iV + i 



By comparing these formulae with that for the quasi- walker, (^), the following simple 
relations are found: 



D. 



bond 



{N + 1 



N + 1 



(115) 



and 



Dhole = D + A 



(116) 



where A is 



M 



k=l 



M 



~^ fc=i 



m. 



(117) 



Eliminating A, the diffusive properties of the hole, the DB and the quasi-walker can be 
related: 



N 

Dhole = {N+ l)Dbond + Jf^^- 



(118) 



These relationships, together with those for the velocities (102 - |106| ), summarize and 
confirm our intuitive picture of the asymptotic properties associated with not only the 
physical objects, the impurity (DB) and the hole, but also the quasi-walker. In the next 
section, we will analyze some special cases in detail, in order to illustrate the use of these 
results. 



4 Special cases 

In this Section we consider some simple, yet interesting, cases. The explicit expressions for 
the velocities and diffusion constants will be transparent, painting a clear picture of the 
physical situation. We will restrict ourselves to two main limits: i) the pure asymmetric 
and symmetric cases and ii) the random symmetric case. In the former, except for jumps 
across the defect bond, the particle-hole exchange rates to the left (right) are uniform, 
given by {W^ )■ Thus, the term "pure asymmetric case" was coined, while the 
"pure symmetric case" consists of a further reduction: = W^. Across the DB, the 
exchange rates (of any particle with the hole) are kept general: q and g', though we also 
consider the limit q = q' Q. For case ii), the exchange rate of any particular particle 
with the hole is randomly chosen, but indepedent of the jump direction - thus the term 
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"random symmetric". The random asymmetric case is extremely interesting, being the 
most general. However, its treatment is quite involved and is beyond the scope of this 
paper. 

As we saw in the preceeding Section, the velocity and diffusion constant of the quasi- 
walker determines almost all interesting quantities of the system. Therefore, we first give 
these expressions in the particular cases. 



4.1 Pure asymmetric and symmetric cases 

To be specific, for the asymmetric case, we have, for all n: 

Wn+l,n = W^ and Wn-l,n = W^ . (119) 

Since all the hopping rates in each direction are identical, the particles in the system are 
no longer distinguishable. It is meaningless to specify "the location of the DB with respect 
to the chain." We can just as easily regard the hole as a particle, suffering biased diffusion 
(if ^ W^) everywhere except across a specific bond (where the bias is q — q'). The 
model reduces to a simple driven system with a single blockage so that its properties 
can be relatively easily understood. 
Define 

s = W^/W^ , (120) 
so that X = (from (|4^)). Using the expressions from Sect. 2.2-3, we then find: 

TT ( As"™' r 1 1 

me[0,N] (121) 



and 



u 



:i - xs) 



m G [0, - 1] 



(122) 



U 



1 - Xs 



m = N 



q — q'X 1 — s 

where U is defined by (43). The normalization constant now takes the form: 



(123) 



UN 



iV + l + 



W^+q' 
q-q'X 



1 - As" 
1 - s 



(124) 



25 



Note that, due to the translational invariance of the pure case, the steady state prob- 
abihty distribution X^^ is independent of b, the DB position. Using (^), we find the 
quasi-walker's velocity V to be: 



V = {W^ - 



1 + 



(125) 



Since the velocities of the physical objects are directly related to this ^ , it is of interest 
to point out two appealing aspects of this expression. The effect of the impurity clearly 
affects V at the level of 0{1/N). Further, it can be completely eliminated by choosing 
q' = and q = W^, where V reduces to the expected result: — W^. In this limit, 
the probability distribution in fact becomes uniform, i.e., X!^ = 1/M. 

A similar expression of the diffusion constant D can be given. It is quite cumbersome 
and not very transparent. So, instead of providing the formula here, let us direct the reader 
to the general case discussed in the last two Sections. However, it is easy to check that, 
if we "eliminate" the DB as above, we recover the expected results: D = {W^ + W^)/2, 
A = 0, Dhole = D, and Dhond = Dhoie/{N + 1)^- The first of these was obtained by 
Derrida [||. The last equation implies that a referrence point (tagged "defect bond") on 
the checkerboard diffuses with respect to the string of particles at a rate N~'^ times slower 
than the hole itself. In other words, the string of particles, as a whole, has a diffusion 
constant given by Dhond- 

The case pointed out in Ref. applicable to the reptation model of a polymer with 
a single drive-insensitive impurity, is given hy q' = q = B. From ( |125| ) we see that the 
impurity reduces the drift velocity by 0{1/N). For this case, the interesting generalization 
of our model lies in the study of chains with a fixed fraction of holes and/or impurities. 

From the "pure asymmetric case", we can simplify further to the symmetric limit, in 
which = = W with arbitrary q and q'. The corresponding physical picture is a 
hole (or a particle) performing simple random walk everywhere except across one bond, 
where it receives an extra "kick" of strength q — q' to the right. Thus, this case can be 
mapped into simple diffusion with open boundaries, coupled to chemical potentials. Solv- 
ing the ordinary diffusion equation with unequal densities at the two boundaries is trivial: 
the result is a linear profile. Indeed, this is reproduced here. Setting s = A = 1 , we find 
the probability distribution: 

6 _ 1 q{N -m) + q'm + W 

- M W+W ' ^^^^^ 

where q = {q + q')/'^ is the average jump rate through the defect bond. The velocity, 
which is related to the current, also takes a simple form: 

V= ""^'^ . (127) 
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However the expression for the diffusion constants remain rather involved. They only 
simplify with the further restriction: q = q' (where clearly the probability distribution is 
uniform and ^ = 0). Then 



D 



1 + N 



(128) 



with A = 0, Dhole = D, and Dbond = D/{N + l)^. 
4.2 Random symmetric case 

In this part we study the case where the W jump rates are symmetric: 

Wn,j = Wj^n 



(129) 



but otherwise randomly chosen out of some distribution. The rates q and q' are also kept 
in general. 

An immediate consequence is A = 1. Using the expressions (44-^) for and and 
introducing 



N 



1 



for convenience, we arrive at: 
and 



1 



/\ Ti I / N-m 

^ 1 + q z 



_j_ 

q' j^^i Wb+m+i,b+m+i-l 



N + 1 
q-q' 



N- 



(130) 



me[0,N] (131) 



(132) 



From the normalization constant 



q{m + 1) + q'{N - m) 



(133) 



M 



^~q 



l + qz 
q-q' 



(134) 
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we obtain the velocity: 



V = ^ . (135) 
l + qz 



This formulae also displays some intuitively understandable properties, e.g., proportion- 
ality to the drive {q — q') and being of 0{1/N) (through z ~ 0{N)). As in the previous 
cases, explicit expressions for the diffusion constants are not particularly illuminating and 
will not be given here. However, in the totally symmetric case, i.e., q = q' , there is simpli- 
fication: 

D = q^±^ (136) 

1 + zq 

with A = 0, Dhole = D, and D^ond = D/{N + 1)^. 

If the ly's are drawn from the distribution piyV), then the expressions for the veloc- 
ity and diffusion constants are valid in the thermodynamic limit N ^ oo with z/N = 
I p( V^)i y-- Finally, we point out the obvious: for a 5 distribution, ( 135 , |136| ) reduce to 
( F27| , |l28|) . 



5 Summary and outlook 

In this paper, we study a generalized asymmetric exclusion process (ASEP) Q to include 
a directional impurity. The simplest realization consists of a random walker (a hole) on 
a one dimensional periodic lattice which is otherwise filled with particles. The rates at 
which the hole exchanges with any particular particle i {i = 1,...,N) is arbitrary and 
direction-dependent (but time-indepedent). In otherwords, the rate of a particle hop to 
the right, Wj+i^j, can be different from those of a hop to the left, Wi-n. However, across 
one special bond (which is fixed in respect to the lattice and not the string of particles), 
which we call the "defect bond", the rates of particle- hole exchange is fixed at q and q' , 
regardless of which particle is involved. In general, the q 's are of course unrelated to the 
Ws. In this sense, there is an "impurity" in the lattice, associated with the mobility of 
the walker, which breaks translational invariance. 

Though the problem appears to have two degrees of freedom - the positions of the 
hole and the defect bond - we show that the associated configuration space is, in fact, 
one-dimensional. There is a subtle coupling of these two objects into a new entity which 
we label the quasi-walker. To be precise, the location of the bond, 6, and the bond-hole 
separation, m, can be combined into a single variable k = 1 -|- (6 — 1)(A^ + \) running 
from 1 to M = A^(A^ -|- 1). A Master equation for Pk{t), the probability to find the 
system in the state k at time t, can be easily written. Being a one-dimensional problem, 
it can be solved and the steady state distribution found explicitly. Further, generalizing 
a method used in the theory of random walks on lattice rings 0,0, (replication) we are 
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able to find implicitly time-dependent information in the steady state, namely, asymptotic 
velocities and diffusion constants associated with both the hole and the defect bond. The 
final expressions, though explicit functions of the g's and W^s, are rather complicated, so 
that it is difficult to gain much insight into their physical content. Thus, we considered 
a few special cases, for which these explicit formulae display intuitively comprehensible 
properties of the system, such as the order of magnitude of the effects of the impurity. 

Turning our focus to physical realizations of this model, we rely on the existing con- 
nection between an impurity-free ASEP and the Duke-Rubinstein reptation model for 
gel-electrophoresis There, a polymer is thought of as a chain of segments, confined 
to a tube of cells. Each cell may contain one or more segments. The links between seg- 
ments are mapped into particles and/or holes in the ASEP in that a link bridging two 
cells corresponds to a particle encoded with "charge" ±1, while a link lying within one 
cell corresponds to a hole. In an external field, the segments move preferentially in the 
direction of the field, a behavior which translates into the "pure asymmetric case" of the 
ASEP. Of course, in general, a polymer will have more than one link lying within a cell, 
motivating the generalization of Derrida's work to a multi-hole ASEP [^]. A more fun- 
damental difference between a periodic ASEP and gel-electrophoresis lies in that open 
boundary conditions are more appropriate for the latter, leading to many changes in the 
string of "particles" in the ASEP. Significantly, a recent work |11] proved rigorously that, 
to leading order in N, the diffusion constant of an open chain coincides with that of the 
periodic chain. Thus, the analysis of a periodic chain [0, 10 1 should be taken more seriously 
than an academic exercise. In a previous publication we have shown that a segment of 
the polymer which is insensitive to the external drive (an impurity) can be mapped into a 
defect bond in the ASEP. Based on the success of the impurity-free ASEP as a model for 
gel-electrophoresis, the work here should be further generalized to include multi-hole and 
multi-impurity cases. The results may be of importance to experiments on the drifting of 
polymers made from a mixture of two different types of monomers. 

Other applications of asymmetric exclusion processes include, e.g., driven diffusive 
systems with two species |jl^ and models of traffic flow [^]. The interesting phenomena 
of phase transitions observed in these systems are, however, present only in two or more 
dimensions. It clearly behooves us to examine our model in dimensions higher than one. 
Another obvious extension of our study is the effect of open boundary conditions. The 
impurity-free ASEP, even in one-dimension, is known to display extremely interesting 
collective behavior such as spontaneous symmetry breaking There is little doubt that 
our study should be regarded as only a small step towards the exploration of the vast 
unknown of non-equilibrium statistical mechanics. 
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